38.2 SCOPE vs STOP

Fit LME models for Prolactin, HPL and GH2.
Fit GEE models for Insulin and HOMA-IR.

  • Outcomes: Prolactin, hPL, GH2, Insulin, HOMA-IR, Glucose (all outcomes log-transformed)
  • Independent variables: Study, GA at sampling
  • Model fitted: LME or GEE model
          log(hormones) ~ Study + GA at sampling + Study:GA at sampling
          LME only: Random intercepts for Study ID and Plate
          GEE only: Independence correlation structure

38.2.1 Global tests

Type III ANOVA on main and interaction effects, and LM model summaries.

PRL (ng/ml)
Chisq Df Pr(>Chisq)
(Intercept) 2.473411 1 0.1 ns
Study 1.757843 1 0.2 ns
GA_samp 8.736573 1 0.003 **
Study:GA_samp 1.550230 1 0.2 ns
hPL (μg/ml)
Chisq Df Pr(>Chisq)
(Intercept) 11.118300 1 0.0009 ***
Study 5.570424 1 0.02 *
GA_samp 33.105635 1 <0.0001 ***
Study:GA_samp 7.350297 1 0.007 **
GH2 (ng/ml)
Chisq Df Pr(>Chisq)
(Intercept) 1.122593 1 0.3 ns
Study 5.234170 1 0.02 *
GA_samp 8.992223 1 0.003 **
Study:GA_samp 5.263896 1 0.02 *
Insulin (uU/ml)
model term df1 df2 F.ratio Chisq p.value
1 Study 1 Inf 0.395 0.395 0.5 ns
3 GA_samp 1 Inf 1.567 1.567 0.2 ns
2 Study:GA_samp 1 Inf 0.731 0.731 0.4 ns
HOMA-IR
model term df1 df2 F.ratio Chisq p.value
1 Study 1 Inf 0.298 0.298 0.6 ns
3 GA_samp 1 Inf 1.230 1.230 0.3 ns
2 Study:GA_samp 1 Inf 0.765 0.765 0.4 ns
Glucose (mmol/L)
model term df1 df2 F.ratio Chisq p.value
1 Study 1 Inf 0.006 0.006 0.9 ns
3 GA_samp 1 Inf 0.003 0.003 1.0 ns
2 Study:GA_samp 1 Inf 0.086 0.086 0.8 ns

38.2.2 Model summary

Ratio of Geometric means (95% CI)

PRL (ng/ml)
Var Effect P
(Intercept) 1.64 (-0.41, 3.67) 0.12
StudySTOP 1.46 (-0.69, 3.61) 0.19
GA_samp 0.20 (0.07, 0.33) 0.004
StudySTOP:GA_samp -0.09 (-0.24, 0.05) 0.21
hPL (μg/ml)
Var Effect P
(Intercept) 1.82 (0.67, 2.96) 0.001
StudySTOP -1.38 (-2.58, -0.17) 0.02
GA_samp 0.21 (0.13, 0.28) <0.0001
StudySTOP:GA_samp 0.11 (0.03, 0.19) 0.008
GH2 (ng/ml)
Var Effect P
(Intercept) -0.69 (-2.09, 0.78) 0.30
StudySTOP -1.61 (-3.17, -0.09) 0.03
GA_samp 0.13 (0.03, 0.22) 0.006
StudySTOP:GA_samp 0.11 (0.01, 0.21) 0.03
Insulin (uU/ml)
Var Effect P
(Intercept) 1.37 (-2.34, 5.08) 0.47
StudySTOP 1.23 (-2.61, 5.07) 0.53
GA_samp 0.14 (-0.10, 0.38) 0.26
StudySTOP:GA_samp -0.11 (-0.37, 0.14) 0.39
HOMA-IR
Var Effect P
(Intercept) 0.11 (-3.61, 3.83) 0.95
StudySTOP 1.08 (-2.79, 4.95) 0.58
GA_samp 0.13 (-0.11, 0.37) 0.29
StudySTOP:GA_samp -0.12 (-0.38, 0.14) 0.38
Glucose (mmol/L)
Var Effect P
(Intercept) 1.71 (0.58, 2.85) 0.003
StudySTOP -0.04 (-1.19, 1.10) 0.94
GA_samp 0.00 (-0.07, 0.08) 0.90
StudySTOP:GA_samp -0.01 (-0.09, 0.07) 0.77

38.2.3 Diagnostics

Model diagnositcs.

Diagnositcs for PRL (ng/ml)

Diagnositcs for hPL (μg/ml)

Diagnositcs for GH2 (ng/ml)

Diagnositcs for Insulin (uU/ml)

Diagnositcs for HOMA-IR

Diagnositcs for Glucose (mmol/L)

38.2.4 Estimated means

Estimated marginal means (or adjusted means) adjusted for GA at sampling.

PRL (ng/ml)
Study GA_samp response SE df lower.CL upper.CL
SCOPE 10 37.84059 14.97325 183.71665 17.33445 82.60489
STOP 10 64.25777 10.79816 25.61462 45.47803 90.79244
SCOPE 16 125.40320 24.78532 46.73448 84.25636 186.64422
STOP 16 122.20695 26.05851 61.88348 79.79445 187.16261
hPL (μg/ml)
Study GA_samp response SE df lower.CL upper.CL
SCOPE 10 48.12118 9.400344 116.3177 32.68197 70.85399
STOP 10 35.67218 1.710800 43.9918 32.38566 39.29222
SCOPE 16 165.09165 11.552861 112.4283 143.71802 189.64396
STOP 16 234.19079 20.961034 179.3138 196.27551 279.43032
GH2 (ng/ml)
Study GA_samp response SE df lower.CL upper.CL
SCOPE 10 1.764828 0.4232374 27.87312 1.0797268 2.884636
STOP 10 1.052072 0.0640941 50.34783 0.9309196 1.188992
SCOPE 16 3.744028 0.3480694 126.03402 3.1148611 4.500281
STOP 16 4.304273 0.5007060 185.72738 3.4216229 5.414615
Insulin (uU/ml)
Study GA_samp response SE df asymp.LCL asymp.UCL
SCOPE 10 15.42851 10.443696 Inf 4.09391 58.14466
STOP 10 17.37172 1.623285 Inf 14.46449 20.86328
SCOPE 16 35.06535 5.621541 Inf 25.61042 48.01086
STOP 16 20.27501 4.212605 Inf 13.49286 30.46618
HOMA-IR
Study GA_samp response SE df asymp.LCL asymp.UCL
SCOPE 10 4.151566 2.7843932 Inf 1.115126 15.456106
STOP 10 3.826518 0.3878844 Inf 3.137036 4.667538
SCOPE 16 9.138816 1.7090499 Inf 6.334418 13.184787
STOP 16 4.200316 0.9587402 Inf 2.685292 6.570106
Glucose (mmol/L)
Study GA_samp response SE df asymp.LCL asymp.UCL
SCOPE 10 5.833361 1.0212224 Inf 4.139067 8.221199
STOP 10 4.947693 0.0651983 Inf 4.821543 5.077144
SCOPE 16 6.006065 0.4569785 Inf 5.173987 6.971957
STOP 16 4.740833 0.1711949 Inf 4.416896 5.088528

38.2.5 Comparisons

Hormone Contrast GA_samp Ratio of Geometric means (95% CI) P-value
PRL (ng/ml) STOP / SCOPE 10 1.70 (0.81, 3.58) 0.2
PRL (ng/ml) STOP / SCOPE 16 0.97 (0.67, 1.42) 0.9
hPL (μg/ml) STOP / SCOPE 10 0.74 (0.50, 1.10) 0.1
hPL (μg/ml) STOP / SCOPE 16 1.42 (1.15, 1.75) 0.001
GH2 (ng/ml) STOP / SCOPE 10 0.60 (0.36, 0.99) 0.05
GH2 (ng/ml) STOP / SCOPE 16 1.15 (0.87, 1.51) 0.3
Insulin (uU/ml) STOP / SCOPE 10 1.13 (0.30, 4.30) 0.9
Insulin (uU/ml) STOP / SCOPE 16 0.58 (0.35, 0.97) 0.04
HOMA-IR STOP / SCOPE 10 0.92 (0.24, 3.48) 0.9
HOMA-IR STOP / SCOPE 16 0.46 (0.26, 0.82) 0.008
Glucose (mmol/L) STOP / SCOPE 10 0.85 (0.60, 1.20) 0.3
Glucose (mmol/L) STOP / SCOPE 16 0.79 (0.67, 0.93) 0.005

38.2.6 Barplots by Study

Barplots in this section may be misleading due to involvement of Study by GA interaction.

38.2.7 Plots for GA